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The Isolated Horizons (IH) formalism, together with a simple phenomenological model for colored 
black holes has been used to predict non-trivial formulae that relate the ADM mass of the solitons 
and hairy Black Holes of Gravity-Matter system on the one hand, and several horizon properties of 
the black holes in the other. In this article, the IH formalism is tested numerically for spherically 
symmetric solutions to an Einstein-Higgs system where hairy black holes were recently found to 
exist. It is shown that the mass formulae still hold and that, by appropriately extending the current 
model, one can account for the behavior of the horizon properties of these new solutions. An 
empirical formula that approximates the ADM mass of hairy solutions is put forward, and some of 
its properties are analyzed. 
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O . I. INTRODUCTION 

U 

In recent years, the introduction of the Isolated Horizon (IH) formalism [J. 0.Q has proved to be useful to gain insight 
into the static sector of theories admitting "hair" 0, [f| . Firstly, it has been found that the Horizon Mass of the black 
hole (BH) , a notion constructed out of purely quasi- local quantities, is related in a simple way to the ADM masses of 
both the colored black hole and the solitons of the theory p|. Second, a simple model for colored black holes as bound 
states of regular black holes and solitons has allowed to provide heuristic explanations for the behavior of horizon 
quantities of those black holes Third, the formalism is appropriate for the formulation of uniqueness conjectures 
for the existence of unique stationary solutions in terms of horizon "charges" |4j. Finally, the combination of the 
Mass formula, together with the fact that in theories such as Einstein- Yang-Mills-Higgs (EYMH) different "branches" 
of static solutions merge, has allowed to have a formula for the difference of soliton masses in terms of black hole 
quantities ■ Many of these predictions have been confirmed in more general situations and for other matter 
couplings UU- For a recent review on IH (including hair) see ^} , and for a review of hairy black holes see • 

In this article, we explore further the consequences of the IH formalism in the static sector of the theory. In 
particular, we explore the behavior of a recently found family of hairy static spherically symmetric (SSS) solutions 
to the Einstein-Higgs system where the scalar potential is allowed to be negative and therefore, the existing 
no-hair theorems [l3| do not apply. In the standard treatment of stationary black holes with killing horizons, one is 
always restoring to several concepts that use asymptotic information very strongly |14| . On the other hand, the IH 
formalism only uses quasi-local information defined on the horizon, allowing it to prove very general results involving 
only these quasi-local quantities. The IH formalism has proved to be generalizable, in the scalar sector, even to the 
non-minimal coupling regime, where the energy conditions required for the consistency of the formalism are much 
weaker [l5| . In the present paper we shall restrict our attention to the minimally coupled case, and for a particular 
form of the scalar potential for which static solutions are known to exist ^3 ■ We will study the one parameter family 
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of solutions (that could be labelled by its geometric radius r&) and compare its properties with those of hairy black 
holes in other theories, such as EYM, where the phenomenological predictions of the IH formalism have been shown 
to work very well |5|,|8j. As we will show, we find that the mass formulae relating BH and soliton ADM masses works 
also well, but the model of a hairy black hole as a bound state of a soliton and a bare black holes exhibits some new 
unexpected features. As we shall see, one needs to slightly modify the model from its original formulation in Ref.0. 
Once this modification is made the model can again explain all the qualitative behavior of the hairy BH solutions. 

The structure of the paper is as follows: In Sec. [H] we review the consequences of the IH formalism for hairy 
solitons and BH solutions. In Sec. IIIII we review the SSS found recently in the Einstein-Higgs sector. Section Hvl is 
the main section of the paper. In it, we show the numerical evidence for the mass formulae and the phenomenological 
predictions of the model. Unlike the EYM case, there are some unexpected features, such as the binding energy 
becoming positive. We then propose a modification of the formalism to deal with such situations. We show that with 
these modifications, the model can still account for the geometrical phenomena found in several theories. In Sec. El 
we explore the situation of the collapse of a hairy black hole and use the model to put bound on the total possible 
energy to be radiated. These results should be of some relevance to full dynamical numerical evolutions of such black 
holes. In Sec. IVII wc propose an empirical formula for the horizon and ADM masses of scalar hairy black holes that 
can also be applied to the EYM case. Finally, we end with a discussion in Sec. IVIII 



In recent years, a new framework tailored to consider situations in which the black hole is in equilibrium ("nothing 
falls in"), but which allows for the exterior region to be dynamical, has been developed. This Isolated Horizons (IH) 
formalism is now in the position of serving as starting point for several applications. Notably, for the extraction 
of physical quantities in numerical relativity and also for quantum entropy calculations P, The basic idea 

is to consider space-times with an interior boundary (to represent the horizon), satisfying quasi- local boundary 
conditions ensuring that the horizon remains 'isolated'. Although the boundary conditions are motivated by geometric 
considerations, they lead to a well defined action principle and Hamiltonian framework. Furthermore, the boundary 
conditions imply that certain 'quasi-local charges', defined at the horizon, remain constant 'in time', and can thus be 
regarded as the analogous of the global charges defined at infinity in the asymptotically flat context. The isolated 
horizons Hamiltonian framework allows to define the notion of Horizon Mass Ma, as a function of the 'horizon charges' 
(hereafter, the subscript "A" stands for a quantity at the horizon). 

In the Einstein-Maxwell and Einstein-Maxwell-Dilaton systems considered originally Q , the horizon mass satisfies 
a Smarr-type formula and a generalized first law in terms of quantities defined exclusively at the horizon (i.e. without 
any reference to infinity) . The introduction of non- linear matter fields like the Yang-Mills field has brought unexpected 
subtleties to the formalism £|. However, one still is in the position of defining a Horizon Mass, and furthermore, this 
Horizon Mass satisfies a first law. 

An isolated horizon is a non-expanding null surface generated by a (null) vector field l a . The IH boundary conditions 
imply that the acceleration k of l a (l a V a l b = K,l b ) is constant on the horizon A. However, the precise value it takes on 
each point of phase space (PS) is not determined a-priori. On the other hand, it is known that for each vector field t„ 
on space-time, the induced vector field X to on phase space is Hamiltonian if and only if there exists a function E to such 
that 5E to = 0(<5, Xt a ), for any vector field S on PS. This condition can be re-written as, SE to — 5a& -(-work terms. 
Thus, the first law arises as a necessary and sufficient condition for the consistency of the Hamiltonian formulation. 
Thus, the allowed vector fields t a will be those for which the first law holds. Note that there are as many 'first 
laws' as allowed vector fields l a =t a on the horizon. However, one would like to have a Physical First Law, where 
the Hamiltonian E to be identified with the 'physical mass' Ma of the horizon. This amounts to finding the 'right k'. 
This 'normalization problem' can be easily overcome in the EM system Q. In this case, one chooses the function 
k = k(cia,Qa) as the corresponding function for the static solution with charges (oa,Qa)- However, for the EYM 
system, this procedure is not as straightforward. A consistent viewpoint is to abandon the notion of a globally defined 
horizon mass on Phase Space, and to define, for each value of n = n (which labels different branches of the solutions), 
a canonical normalization that yields the Horizon Mass M^™ ' for the n a branch The horizon mass takes 

the form (from now on we shall omit the no label), 



II. CONSEQUENCES OF THE ISOLATED HORIZONS FORMALISM 
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with ta the horizon radius. Here fi(r&) is related to the surface gravity as follows P(ta) — 2rAK(r&). 

Furthermore, one can relate the horizon mass Ma to the ADM mass of static black holes. Recall first that general 
Hamiltonian considerations imply that the total Hamiltonian, consisting of a term at infinity, the ADM mass, and a 
term at the horizon, the Horizon Mass, is constant on every connected component of static solutions (provided the 
evolution vector field agrees with the static Killing field everywhere on this connected component) {HQ In the 
Einstein- Yang-Mills case, since the Hamiltonian is constant on any branch, we can evaluate it at the solution with 
zero horizon area. This is just the soliton, for which the horizon area a a, and the horizon mass Ma vanish. Hence 
we have that = M so \. Thus, we conclude 

M so i = Madm - M A , (2) 

Thus, the ADM mass contains two contributions, one attributed to the black hole horizon and the other to the outside 
'hair', captured by the 'solitonic residue'. The formula (0), together with some energetic considerations 0, lead to 
the model of a colored black hole as a bound state of an ordinary, 'bare', black hole and a 'solitonic residue', where 
the ADM mass of the colored black hole of radius r& is given by the ADM mass of the soliton plus the horizon mass 
of the 'bare' black hole plus the binding energy: 

Madm = M so] + M A = M A + M sol + E hind , (3) 

with -Ebind = Ma — M A . Simple considerations about the behavior of the ADM masses of the colored black holes and 
the solitons, together with some expectations of this model (such as demanding for a non-positive binding energy) 
give raise to several predictions about the behavior of the horizon parameters Among the predictions, we have: 

i) The absolute value of the binding energy decreases as ^a increases. 

ii) /3(rA), as a function of r&, is a positive function, bounded above by /3(o)(^a) = 1- 

iii) The curve /3(r), as functions of r intersect the r = axis at distinct points between and 1, and never intersect. 
Finally, 

iv) The curve for /3, for large value of its argument, becomes asymptotically tangential to the curve /3(o)(rA) = 1- 
One of the features of these solutions in, say, EYM is that there is no limit for the size of the black hole. That 

is, if we plot the ADM mass of the BH as function of the radius r& we get an infinite number of curves, each of the 
intersecting the rA = line at the value of the soliton mass, and never intersecting each other. 

The purpose of this paper is to test the mass formula J5J, for the scalar hairy solutions found in Ref. 01 an d also 
to confront the predictions i)-iv) (obtained for the colored EYM BH model p|), with the corresponding properties for 
the scalar hairy BH. In the next section we will review the scalar hairy solutions, and afterwards we shall study their 
horizon properties. 



III. SCALAR SOLITONS AND BLACK HOLES IN EINSTEIN-HIGGS THEORY 



Let us consider the theory of a scalar field minimally coupled to gravity described by the total action: 



Stotb^, 4>] 



R 

16^ 



-(v^)(v^) + y(0) 



d 4 x 



(units where Go = c = 1 are employed). The field equations following from the variation of the action Q are, 



G 



and, 
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dV{(j>) 
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It is well known that asymptotically flat static spherically symmetric solutions representing black holes solutions 
to the Einstein-Higgs equations do not exist if the scalar matter satisfies the weak energy condition (WEC) due 
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to the existence of the so called scalar no- hair theorems jT3|. Recently, numerical evidence for the existence of 
asymptotically flat and static spherically symmetric solutions representing scalar hairy black holes (SHBH) and scalar 
solitons (scalarons; hereafter SS) have been found in theories represented by the action Q and with a scalar potential 
non-positive-semidefinite ^| given by the asymmetric potential, 



v(<f>) 



4(771 + J72) 



2771772 



(7) 



where A, T)i and a are constants. For this class of potential one can see that, for 771 > 2ij2 > 0, <fi = a corresponds to 
the local minimum, 4> = a + r\\ is a global minimum and (f> = a + 772 is a local maximum (see Fig. 1). The key point in 
the shape of the potential, V(6), for the asymptotically flat solutions to exist, is that the local minimum Vj]^ 1 = V(a) 
is also a zero of V(<f>) (to see [l2j for an detailed analysis for the existence of these solutions). Moreover, V{4>) is not 
positive definite (we assume A > 0), which leads to a violation of the WEC and therefore the scalar no- hair theorems 
[l3| can not be applied to this case. 




FIG. 1: 

Qualitative shape of the scalar-field potential V(<f>) as given by Eq. £J used to construct the asymptotically flat black hole 

and soliton solutions. 

In order to describe the asymptotically flat SHBH and SS, we use a standard parametrization for the metric and 
the scalar field describing spherically symmetric and static spacetimes 



ds 2 = -11- 



2m(r) 



dr 2 +r 2 dfl 2 , 
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For SHBH we demand regularity on the event horizon r A which implies the conditions, 



(8) 
(9) 

(10) 
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(d r m) A = ^r 2 A V A . 



(1 - 87rr A y A ) 

For SS we impose regularity at the origin of coordinates r = 0, 

m(0) = 0, 8(0)= J , 0(0) -0o, (cW) o = 0. 



(11) 



(121 



where So and 0o are to be found such as to obtain the desired asymptotic conditions. In addition to the regularity 
conditions, we impose asymptotically flat conditions on the spacetime for SHBH and SS: 



m(oo) = M A dm, S(oo) = 0, 0(oo) = C 



(13) 



■5 



Above, the value tfioo corresponds to the local minimum of V(<p). A/adm is the ADM mass associated with a SHBH 
or SS configuration. For a given theory, the family of SHBH configurations is parametrized by the free parameter 
A A which specifies the area of the black hole horizon. Therefore for SHBH, Madm = ^adm(^a), or equivalently 
Madm = AfADM( r A) since in our coordinates the horizon area A A = 47rr A . The value 4> A = 4> A ( r A ) is a shooting 
parameter rather than an arbitrary boundary value which is determined so that the asymptotic flat conditions are 
satisfied. On the other hand, for SS the value cj>o is the shooting parameter, and the corresponding configuration is 
characterized by a unique M A ^ M . 

Finally, the surface gravity of a spherically symmetric static black hole can be calculated from the general expression 
of spacetimes admitting a Killing horizon 



lly s 

m 



1/2 



(14) 



where V 2 stands for the Laplacian operator associated with the stationary metric and a = (dt, dt) is the norm of the 
time-like (static) Killing field which is null at the horizon. For the present case, a 
From the above formula, one can obtain the following useful expression 
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For the election of the parametrization of the metric JSJ we have 
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«(r A ) = —e d{r *> [1 - 2{d r m) A ] . 
2r A 

Introducing in (|16|l we obtain the final expression for the surface gravity of the SHBH 

1 



n(r A ) 



2r A 



[1 - 8icr 2 A V A ] 



(15) 
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In the next section, we shall analyze these solutions from the perspective of the IH formalism. 



IV. MASS FORMULAE 



Let us now turn to the straightforward application of the IH formalism mentioned in the Sec. II to the case of SHBH 
and SS in the Einstein-Higgs theory with action given by (0J and V{4>) by Eq. (0) 12]. As in Ref. 01 > we shall 
take the specific values 771 = 0.5, 772 = 0.1 and a — 0; all the quantities (e.g., Madm and r A ) have been rescaled as 
appropriate using l/\/A as a length-unit. 

The first consequence coming from the IH formalism is that the horizon mass associated to the SHBH takes the 
form 1 , 

1 r A 

M A (r A ) = -J /?(r)dr, (18) 

where /3(r A ) — 2r A n{r A ) [the value of n(r A ) is given by (|17[l ]. We have dropped the subindex (n) in the expression 
(118(1 because in the Einstein-Higgs system considered here there is only one branch of static spherically symmetric 
SHBH labelled by its horizon radius r A (the corresponding scalar configurations do not have nodes). Additionally, 
there is another branch of static spherically symmetric BH given by the family of Schwarzschild BH's labelled by its 
corresponding horizon radius r A and with horizon mass 

M A chwarz (r A ) =r A /2 . (19) 



1 We remind the reader our choice of units Go = c = 1 . 
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FIG. 2: 

The ADM mass (solid lines), the horizon BH mass (dash-dotted lines), and the mass of the Schwarzschild solution 
(dashed- lines) plotted as functions of r& (first panel). The second panel depicts similar quantities using logarithmic scales to 
appreciate better their behavior for small ta- The soliton mass M ao i ~ 3.827. 
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FIG. 3: 

P{ta) = 2k rA is plotted as function of ta- Note that it approaches asymptotically the value « 1.24. Here /3(0) ~ 0.324. 

The second consequence coming from the IH formalism is that on the entire branch of SHBH we can expect the 
following identity to be true 

Madm = M so i + M A , (20) 

where M so \ is the ADM mass of the SS obtained taking the limit r& — > of the branch of SHBH and Madm is the 
ADM mass corresponding to the SHBH with horizon radius ta, and the horizon mass is given by Ijl8|) . Thus, in 
a similar way to the EYM theory, the total ADM Mass of the solution contains two contributions, one attributed 
to the horizon of the SHBH and the other to the outside 'hair', captured by the SS. We have performed numerical 
explorations for SHBH for a large range of values of the horizon radius (in normalized units) and have checked the 
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identity l|2UI) . We have found complete agreement within the numerical uncertainties. This can be seen in Fig. |3 
where the identity was checked up to ta = 250. 

Figure depicts the behavior of (i{r&). Unlike the EYM model, where (3 w 1 for large r&, in this model j3 ~ 1.24 
asymptotically. 

Figure 0] shows an example of a BH solution with large ta (for a small BH see fig.2 of Ref. 
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FIG. 4: 

Large-black- hole configuration constructed with V{4>) as given by Eq. with parameters 771 = 0.5, J72 =0.1, a = 0, and 
ta = 150/\/A, 0a ~ 0.26111. The upper panels depict the scalar field and the mass function respectively. The latter 
converges to Madm ~ 93.096/\/A. The lower panels depicts the metric potentials (the first is a zoom of the second): y/—gtt 
(solid line), yjg rr (dashed line), e s (dash-dotted line) and 8 (dotted line). 



A. A physical model of SHBH 



The extrapolation of the model of a hairy black hole as a bound state of an ordinary, 'bare', black hole and a 'solitonic 
residue' (first applied successfully to the colored BH's in the EYM theory) does not apply directly to the SHBH 
because the straightforward generalization of the formula as 

M ADM = M sol + M A = Mi° hwarz + M sol + E hind , (21) 

where -E^ind = Ma — M A chwarz , to our case, has the problem that the binding energy changes sign, becoming positive 
for BH larger that ta ~ 30, and then increasing in absolute value as ta gets larger. That is, 



#bind ~ Br a 



(22) 
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(where B is a constant whose value depends on the specific model) for > 1, which is contrary to the expected 
feature of a negative binding energy as in the EYM case (i.e. the prediction (i) mentioned above will not be satisfied). 

In order to appreciate the origin of the failure of this feature, let us recall that we can write Madm i n terms of the 
Schwarzschild mass and the mass of the 'hair' as, 

M AD M(r A ) = Ml chwarz (r A ) + M hair (r A ) , (23) 

where 

M hail (r A ) = - r 2 T{dr , (24) 

J l-A 

By rescaling the r— coordinate in terms of r A , we can rewrite the mass of the hair. It becomes then, 

/oo 
x 2 f}dx , (25) 

with 

x = — , (26) 
r A 

T* = r|T* , (27) 
2m(r)\ (d r (f) 



(28) 



Now, for x 3> 1, the integral in Eq. l|25|l becomes almost independent of r A , and in fact the numerical analysis 
provides the following value 

M h air(r A >l)~Br A . (29) 

where B « 0.12. Now, since Ai^ chwarz (r A ) = r A /2 we have then 

M ADU {r A > 1) ~ Cr A . (30) 

where now C « 0.62. Therefore we conjecture that C is a constant that depends of the matter-theory involved. 
For the EYM case, one can easily show that the scaling properties of the hair contribution of the energy makes the 
equivalent of the integral of Eq. l|25|) to behave like l/r A rather than r A . Therefore, for the EYM case, C = 1/2. 
As we now show, this subtle difference in both theories makes that the binding energy expression used in the EYM 
cannot be used straightforwardly for the Einstein-Higgs theory analyzed here. 
From and (gSJ we find that 

E hind (r A > 1) - M ha ir(r A > 1) - M sol ~Br A - M soh (31) 

then £"bind scales as r A when r A 3> 1 (remember that M so \ is a constant) [T3 |. It is clear that the sign of the binding 
energy as such defined changes sign and grows with the size of the black hole. 

To deal with this situation we now proceed to propose a modification of the model of a hairy black hole as a bound 
state in order to adapt it to the more general case. Our proposal consists in redefining the binding energy as, 

Kin W d(^A) = E hind (r A ) — B r A = M A (r A ) - Ml chwarz (r A ) - Br A , (32) 

By this procedure we have 'renormalized' the binding energy by subtracting the divergent term. That is, the new 
expression for the binding energy is 

S(rA) = M A {r A ) - Cr A . (33) 

This definition shares now exactly the same properties as the original expression for the EYM. Thus, it vanishes at 
r A = 0, and decreases monotonically to the negative value —M so \. It is interesting to note that this new definition 
reduces to the old one for the EYM case, since as we previously remarked Ceym = 1/2- 
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FIG. 5: 

Both binding energies are plotted. The 'old' binding energy -Ebind^A) is shown to become positive and approaches 
asymptotically a straight line. The 'new' binding energy E^™ d (r/^.) is also plotted (continuous line), showing the expected 

behavior. 



The formula (|21(1 can now be reformulated as, 



Madm(fa) = M so i + M A (r A ) 



(34) 



schwarz 



Ml 
Cr A 



M so i -f 



hM sol 

77 new 



Br A 



Thus, we would get the same structure for the ADM mass of the hairy black hole where now the mass of the 'bare 
black hole' would be equal to M A = (1/2 + B)r A = Cr A . However it is not clear what the origin of this extra term 
(Br A ) is, and we could very well have assigned it to the soliton mass to form a new 'solitonic residue with mass 
M so \ + Br A (whose interpretation however seems somewhat obscure). We must explore more in order to decide which 
interpretation is best suited for our model. As we stressed, the constant C depends on the theory considered; for 
the EYM and EYMH theories, C = 1/2. It would be interesting to explore (in addition to the current analysis, see 
below) whether other theories admitting hair posses values of C different from 1/2. 

To end this section, let us rewrite the form of the hairy mass in terms of the old binding energy, in order to 
understand the behavior of the scalar system. First, let us note using Eqs. (|21(l and l|23|l that 



M hail (r A ) = Madm^a) 



j^j- schwarz 



(r A ) = M so i + E ] 



bind 



so the binding energy is 



-Ebind = MhairfTA) ~ M ao \ = - 



r 2 Tldr 



r 2 T \ dr 



o 



where T Q t is the stress-energy tensor of the solitonic regular solution. This equation can be rewritten as, 



/ r 2 (T/-T * t )dr- / r 2 T \dr 



(35) 



(36) 



(37) 



Here we can identify the first term as the difference between the hair of the BH and the "hair" of the soliton. Of 
course we are comparing the quantities (the integrands) that live on different manifolds, but the total integral is well 
defined. What happens in the EYM case is that both stress tensors behave very much alike, for the exterior region 
(r > r A ) and for large values of r A , and thus the only term that contributes is the second one, that gives the ADM 
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FIG. 6: 

This figure illustrates a physical process where an initial configuration with an isolated horizon Ai n is perturbed and the final 

state contains another isolated horizon Afl n . 

mass of the soliton (recall that r& ^> 1, that is for BH's much larger than the characteristic size of the soliton (of order 
one in this dimension-less units), so the integral captures most of the soliton mass). In the scalar field case, the fact 
that the binding energy is proportional to r&, for large black holes, is captured by the fact that the BH contribution 
to the first term in l|37|l is dominating. It would be interesting to explore this issue in other gravity-matter systems. 

V. INSTABILITY AND FINAL STATE 

The next question we want to consider has to do with the following situation. Consider the case where a hairy black 
hole of geometrical radius r& is slightly perturbed and therefore it decays. The final state will be, one expects, a 
black hole that in its near horizon geometry resembles the Schwarzschild solution, with the scalar field taking the 
value where the potential has a local minima and vanishes. This means that in this process the "scalar charge" at 
the horizon, namely the value (/>a must change. One can make an argument similar to the one in Ref. |a] to conclude 
that, in that situation, the horizon must grow in the process and therefore, the available energy to be radiated can 
not all be radiated to infinity; part of it must fall into the black hole. Let us now recall the estimate for the upper 
bound of the total energy to be radiated. 

The first step is to assume that the process illustrated in Fig. El takes place. Then, we assume that in the initial 
surface there was an isolated horizon Aj n and after the initial unstable configuration has decayed, with part of the 
energy falling through the horizon and the rest radiating away to infinity, we are left with a horizon Afi n of a hairless 
black hole (with r A n > r^). If we denote by Ej+ the energy radiated to future null infinity I + , and given that the 
ADM energy does not change in the process, we have 

Madm - Ma (?"a ) + M sol = Ml chwarz (r A n ) + E x+ , (38) 

which can be rewritten as, 

M sol + E™ nd = Ml chwarz (r A n ) - Mi chwara (r£ ) + E I+ . (39) 

On the right-hand-side note that the first two terms can be identified with AM A , namely the change in (bare) horizon 
mass, while the second term corresponds to the radiated energy. Thus, it is natural to identify the quantity on the 
left as the available energy -E a vaii on the system. We can then write, 

£avaii = M ADM - M£ hwOT (r£) = M sol + E bind (r%) = M sol + Cd(^) + B r% (40) 
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There are several comments regarding this quantity. First, we note that there is a qualitative change in the behavior 
of -Eavai^^A) as function of the initial horizon radius, as in the EYM case. Its functional dependence is very similar to 
the binding energy since they differ only by the soliton mass. In the EYM case the available energy was equal to the 
soliton mass when there was no initial black hole (there is no energy used in binding the BH), and decreases as the 
radius increases. For very large black holes, the available energy goes to zero. For the scalar case under consideration 
here, we have a different behavior. The available energy decreases for small black holes but starts to increase and 
grows linearly with r& ■ The fact that in the EYM case the available energy went to zero for large BHs was interpreted 
as meaning that those black holes were 'less unstable'. This expectation was confirmed by the fact that the frequencies 
of the linear perturbations was decreasing with the radius of the initial BH 0, 0] . It is natural then to ask the same 
question for the scalar black holes. We have computed the frequencies of the (single) unstable mode ip(t, r) = x(r)e t,yt 
present (where a 2 turns to be always negative), as a function of the horizon radius and plotted it in Fig. \7\ (where 
ip(t,r) represents a linear perturbation of <f)(r); see Ref.[T^| for the details). 
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FIG. 7: 

The unstable-mode frequency uj = y/—a 2 of the perturbed BH and soliton (rA = 0) is plotted as a function of the horizon 
radius ta- Note that the frequency decreases as the horizon radius increases. 

As can be seen from the figure the frequencies still decrease as the black holes become larger, which is the same 
behavior observed in the EYM case. It is convenient then to reconsider the meaning of 'less unstable'. In the scalar 
field case considered here, numerical investigations of the dynamical evolution of the soliton as initial state show that 
the system is unstable 18]. The dynamical evolution of the system depends on the sign of the initial perturbation 
on the extrinsic curvature. For one sign of the perturbation, the system collapses and forms a black hole with a final 
isolated horizon,while for the other sign the system expands as a domain wall and gets therefore radiated to infinity 
(for details see [l8j). One should then expect that the dynamical evolution of slightly perturbed hairy black holes 
will show a similar qualitative behavior. In that case, for one sign of the perturbation one might expect the situation 
considered before, namely that the scalar field collapses and the BH grows. For the other sign, one can imagine that 
there could be, in some situations, an expanding wall that radiates away while leaving a "naked black hole" . The 
pressing question in whether, in that case, this residue would be a Schwarzschild like or an AdS like black hole. This 
question arises since, for the soliton collapse in the case of the expanding wall, the region around the origin resembles 
an AdS spacetime with an effective negative cosmological constant generated by the (non-positive) potential. One 
might need in that case a new interpretation of the formalism. 

It is our belief that one needs to clarify what the criteria should be for regarding the system as slightly unstable 
or very unstable, other than the frequency of its perturbations. This and a full clarification of the nature of the 
resulting bare black hole could be achieved whenever full numerical simulations of dynamical evolution staring from 
scalar hairy black holes become available. 

Let us return our discussion to Eq. I|4l)|l . The first thing to note is that due to the characteristic behavior of 
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these solutions, for horizons larger than ta ~ 30, the horizon mass ol the hairy black hole becomes larger than the 
Schwarzschild horizon mass of the same radius. This is also the point at which the binding energy becomes positive. 
One can thus speculate that the black holes of this radius and larger will have more violent collapses with a larger 
fraction of the available energy radiated away. Finally, from Eq. H4l)|) one could interpret that again, the term (-Br a ) 
could be associated to the soliton mass to form a solitonic residue that, together with the new binding energy allows 
us to have the same qualitative features of the heuristic model of [5j. Again, a more detailed analysis will have to 
wait for the numerical investigations of the fully dynamical process. 



VI. AN EMPIRICAL FORMULA 



The non-linear behavior of the Einstein-Matter equations and the non-trivial relations between the masses and the 
horizon radius posses a challenge to obtain an analytical formula for Madm(^a)- One should expect that there is in 
general no closed analytical formula for the masses of hairy BH. 

We have discovered that the following empirical formula reproduces the qualitative and quantitative features of the 
numerical analysis 



M ADM( r A) = 



M H , 



sol 



Ms. 



sol 



D 

2C 



,,2 
A 



D 



Cr A 



A/ S ol 



(41) 



where 



D = 



2C(3 M sol 
2C-/3 



M so , = 



(3 M sol 
2C-[3 



Po being the value of /? at ta = 0. 

The formula (|41(1 has the following nice properties: 



(42) 



1. MX7 M (0) = AL 



'ADM 
rcmp 



sol • 



2. Af AD ^j (rA ) is a monotonically increasing function of ta- 

3. For large r A , M A dm CrA- 

4. The relative error between Aff D ?. and the numerical one is less that 10%. These errors become very small for 
small and large 



5. One can define an empirical binding energy by using Ef- (t"a) 



M emp 
ADM 



M BO \ — CrA [where the first two 

terms provide the empirical horizon mass; here we are using Eq. (|33|l ]. This formula reproduces very well the 
numerical results. 



6. One can then obtain a fit for j3 as follows 



P cm (r A ) = 2 



drA 



2r A C 2 + D 



M 



sol 



§,) A + C 2 r 2 A + Dr A 



(43) 



1 



Cr A 
M Bo] 



F 



This formula reproduces the qualitative shape of the numerical 0, such as its exact value at the origin, its 
monotonically increasing behavior and its asymptotic value 2C for large ta- 



7. For the Schwarzschild case (C — 1/2, M f 



sol 



HrA)EU3(r A )E0 



0), one obtains the expected results: 



MXdm(^a) 



ta/2, 
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Clearly by adding terms of the form (1 < 9 < 2) inside the square root of Ea. (|41|) one could improve the fit 
between the numerical results and the analytical formula. In Figure [3] we compare between the empirical formula and 
the numerical values of the hairy scalar black holes, for the ADM mass, screened surface gravity j3 and the binding 
energy. 

The empirical formula can be used also for the EYM case with C = 1/2 and the corresponding values of M so \ and 
D. Figure [5] compares the numerical values of the EYM n — 1-branch with those obtained from the empirical formula. 

We conjecture that the empirical formula can work also for different n, by using their corresponding values M s "j, 
and Pq . Moreover, we also speculate that such a formula can hold for other theories admitting hair, such as in the 
Einstein-Skyrme and Einstcin-sphalcron models. It remains to be investigated what are the values of C, M so \ and (3q 
for such cases. 

Now, we can further use the first law of thermodynamics SM — kSAa/(Stt) for Af AD j^, and obtain the following 
prediction 

«("H. + 5g) + («> 

Since the properties 1)— 7) show that the analytical results obtained from Eq. 14111 work particularly well for large 
and small r&, the most reliable consequence of (|44[1 is a remarkable simple relation between the surface gravity and 
the ADM mass for sufficiently large hairy black holes: 

K M ABM ^C 2 . (45) 

Note that Eq. (|45|l is consistent for the Schwarzschild case (C = 1/2, D = 0, A/adm = ^a/2), where the identity 
kMadm = 1/4, holds exactly. 
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FIG. 8: 

The three panels depict the ADM-mass, /3, and binding energy E^™ d , respectively, as a function of the horizon radius. The 
solid lines correspond to the values obtained from a numerical analysis and the dashed lines were obtained from the empirical 
formulae described in the main text. Note the good qualitative behavior of the empirical formulae. Remarkably good fits to 
the more precise numerical values are obtained for small and large r&. The values for the empirical formula are C ~ 0.62, 

M so i « 3.827 and /3 ~ 0.324. 



We have performed a non-exhaustive analysis of the solutions with respect to variations of some of the parameters 
of the scalar potential Eq. J7J. Notably, we have computed the effect of the variation of 772 on the global quantities. 
It is to note that changing 772 modifies the potential barrier between the global and the local minimum. In fact, the 
closer the value 772 to 771 /2, the less negative is ^(0 + 771) = \-qf [2t]2 — r/i] /12, and therefore the potential approaches 
the conditions where the no-hair theorems apply. The details of the solutions then depend in a non-trivial fashion 
between the interplay of the negative global minimum (in order to avoid the applicability of the non-hair theorems) 
and the height of the potential barrier. 

Figure I1UI depicts different global quantities as a function of the horizon radius for five different values of 772 ■ The 
soliton mass (r& = 0) as well as the ADM and horizon masses (for large r&) tend to increase with 772 . Remarkably, 
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FIG. 9: 

Same as FiglHl for the Einstein- Yang-Mills theory (n = 1 colored black holes). Here the values for the empirical formulae are 

C = 1/2, Msoi « 0.828 and fa » 0.126. 

the empirical formulae continue to provide reasonable good results by changing the corresponding values of their 
parameters P m :— (M so i, C,/3q) . The quality of the fit to the numerical values can be appreciated by the dashed 
curves of Fig. 1101 which were computed with the empirical formulae. 

VII. DISCUSSION 



Let us first summarize our results. By solving numerically Einstein's equations for static solutions of a self-gravitating 
scalar field, we have analyzed the behavior of several spacetime quantities as functions of the black hole horizon radius. 
We have found that the ADM mass of the spacetimes exhibits two types of behavior: it is similar to other 'hairy' 
theories for small black holes, but its behavior changes dramatically for large black holes. In particular the ADM mass 
of large BH scales not as ta/2 as in other theories (EYM, EYMH, etc), but the proportionality constant (with respect 
to horizon radius) takes a different value depending on the form of the potential (C w 0.63 for rj 2 = 0.1). In this article 
we have analyzed the consequences of this fact for a model based on the isolated horizons formalism. In such a model, 
a hairy black hole is viewed as a bound state of a soliton (which we have) and a 'bare black hole'. The binding energy 
is found to be negative in EYM and EYMH, but in our case, for large BH, the binding energy becomes positive and 
grows linearly with r<\ 2 ■ This fact leads to several possibilities. We have seen that it is possible to modify the original 
model by 'renormalizing' the binding energy in such a way that the newly defined energy has the same qualitative 
behavior as in the EYM system. The price one has to pay is the need to reinterpret either a new 'solitonic residue', or 
a new bare black hole. As a first attempt towards giving a definite answer to this question, we analyzed the frequency 
of the unstable mode of the linearized perturbation, and found that the behavior is the same as in EYM. This suggests 
that the proper physical interpretation is still unclear and that further numerical dynamical investigations are needed 
to fully settle the question. In particular, the two different regimes of the theory might have some consequences in 
the dynamical evolution of slightly perturbed BHs, where one could conjecture a different qualitative behavior for 
small and large black holes, regarding the endpoint of evolution and the nature of the bare black hole to which the 
solution settles. We have also conjectured that the constant that fixed the proportionality between ADM mass and 
horizon radius for large BH's is a theory-dependent constant, which would in particular imply that axi-symmetric 
non-spherical BH solutions to the gravity-scalar field system would have the same asymptotic behavior, for each given 
potential. It would be worth studying other gravity-matter systems, such as non-minimally coupled scalars, to see 
whether they posses a different proportionality constant (work is in progress in these directions). 

We have shown also that a very simple heuristic analytic formula captures the essential qualitative behavior of the 
ADM mass of the hairy scalar BH's, specially for small and large values of the horizon radius. We have conjectured 



2 Recently, another system in 5 dimensions was shown to posses a positive binding energy as well 
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FIG. 10: 

Panels 1-4 depict the ADM- mass, horizon mass, j3 and the binding energy respectively as a function of 7"a- The solid lines 
were obtained from the numerical analysis while the dashed lines were computed using the empirical formulae. The lines are 
associated with the five different values used for 772 = 0.1, 0.11, 0.12, 0.13, 0.14 with 771 = 0.5 fixed. As seen from bottom to top 
(for large ta) the plots of panels 1-3 correspond to 7/2 in increasing order (in panel 4 the order is reversed). The values of the 
parameters P V2 = (M so i, C, /3o) used in the empirical formulae are P0.1 ~ (3.82,0.62,0.32), Po.11 ~ (5.22,0.66,0.24), 
P0.12 « (7.47,0.72,0.17) P0.13 « (11.78,0.82,0.11) P .i4 « (23.51,1.01,0.06). 

that such formula can also be useful for EYM and more general hairy black holes. It remains a theoretical challenge 
to fully understand the origin of such simple formula. 

Perhaps the most important conclusion from the present work is the lesson that hairy black holes for different 
matter systems exhibit new, and sometimes, unexpected behavior. This also point out to the need of a proper and 
deeper understanding of the reason why the heuristic hairy black hole model works so well for the system that it does, 
and whether the phenomenological modifications we have proposed here stand the test of full numerical investigations. 

Acknowledgments 

We would like to thank A. Ashtekar and D. Sudarsky for discussions. This work was in part supported by grants 
D G APA-UN AM IN122002, and INI 19005. U.N. acknowledges partial support from SNI, and Grants No. 4.8 CIC- 
UMSNH, No. PROMEP PTC-61 and No. CONACYT 42949-F. 



[1] A. Ashtekar, C. Beetle, O. Dreyer, S. Fairhurst, B. Krishnan, J. Lewandowski and J. Wisniewski, "Isolated horizons and 
their applications," Phys. Rev. Lett. 85, 3564 (2000), gr-qc/0006006 



16 



[2] A. Ashtekar, C. Beetle and S. Fairhurst, "Mechanics of Isolated Horizons," Class. Quant. Grav. 17, 253 (2000), 

gr-qc/9907068 A. Ashtekar and A. Corichi, "Laws governing isolated horizons: Inclusion of dilaton couplings," Class. 

Quant. Grav. 17, 1317 (2000), gr-qc/9910068 
[3] A. Ashtekar, S. Fairhurst and B. Krishnan, "Isolated horizons: Hamiltonian evolution and the first law," Phys. Rev. D 62, 

104025 (2000), gr-qc/0005083 

[4] A. Corichi and D. Sudarsky, "Mass of colored black holes," Phys. Rev. D 61, 101501(R) (2000), gr-qc/9912032 

A. Corichi, U. Nucamendi and D. Sudarsky, "Einstein- Yang-Mills isolated horizons: Phase space, mechanics, hair and 

conjectures," Phys. Rev. D 62, 044046 (2000), gr-qc/0002078 
[5] A. Ashtekar, A. Corichi and D. Sudarsky, "Hairy black holes, horizon mass and solitons," Class. Quant. Grav. 18, 919 

(2001), gr-qc/0011081 

[6] B. Kleihaus and J. Kunz, "Non-Abelian black holes with magnetic dipole hair," Phys. Lett. B 494, 130 (2000), 
hep-th/0008034 

[7] A. Corichi, U. Nucam endi and D. Sudarsky, "A mass formula for EYM solitons," Phys. Rev. D 64, 107501 (2001), 
|arXiv : gr-qc /0 106084 

[8] B. Kleihaus, J. Kunz, A. Sood and M. Wirschins, "Horizon properties of Einstein- Yang-Mills black hole," Phys. Rev. D 

65, 061502(R) (2002), arXiv:gr-qc/0110084 
[9] R. Ibadov, B. Kleihaus, J. Kunz and M. Wirschins, "New black hole solutions with axial symmetry in Einstein- Yang-Mills 

theory," Phys. Lett. B 627, 180 (2005) arXiv : gr-qc/0507110 
[10] A. Ashtekar and B. Krishnan, "Isolated and dynamical horizons and their applications," Living Rev. Rel. 7, 10 (2004), 

|arXiv : gr-qc/ 0407042 

[11] M. S. Volkov and D. V. Gal'tsov, "Gravitating non-Abelian solitons and black holes with Yang-Mills fields," Phys. Rept. 

319, 1 (1999), |arXT77he p-th/9810070 
[12] U. Nucamendi and M. Salgado, "Scalar hairy black holes and solitons in asymptotically flat spacetimes," Phys. Rev. D 

68, 044026 (2003), |arXiv7gr -qc/030 1062"1 
[13] M. Heusler, J. Math. Phys. 33, 3497 (1992); D. Sudarsky, Class. Quantum Grav. 12, 579 (1995); J.D. Bekenstein, Phys. 

Rev. D51, R6608 (1995). 
[14] M. Heusler, Black Hole Uniqueness Theorems, Cambridge Univ. Press, Cambridge (1996) 

[15] A. Ashtekar, A. Corichi and D. Sudarsky, "Non-minimally coupled scalar fields and isolated horizons," Class. Quant. Grav. 

20, 3413 (2003), arXiv:gr-qc/0305044 
[16] Another example of hairy BH's with the property that the masses Mhair scale also linearly with ta are the hairy BH's of 

the Einstein- Skyrme system. To see this we take the equations of motion for the Einstein-Skyrme model as described in 

the subsection (7.4) from the reference |Tl|: specifically, we take the equation (7.37): 



dm(r) 
dr 



N 



+ sin x 



2 sin 2 y 
r H i 



sin 2 x 



(46) 



where the metric for static and spherically symmetric configurations is 
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with 



N(r) =1- 2m(r)/r, 
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and the Skyrmion field x( r ) depending on the dimensionless coordinate r; to continuation we define a new coordinate as 
x — r/rA and after integration rewrite l|460 as 
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The behavior for x( x ) can see from left panel in the fig. 14. of such reference. From equation (7.41) we see the asymptotic 
behavior for x( x ) as x ~~ * 00 (with a a constant): x ~ ax" 2 . Altough this analysis shows that Madm(?"a) grows linearly 
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to configurations with horizon radius ta < r^^i^) (here n = 4-kGo/ 2 is the coupling constant of the theory and r~A a:L {K) 
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